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Abstract 

Generally, the measurement process consists in coupling a system to a detector that can give a 
continuous output. However, it may be interesting to use as a detector a system with a discrete 
spectrum, especially in view of applications to quantum information. Here, we study 1) a two- 
level system measuring another two- level system (qubit); 2) a generic system measuring a qubit; 
3) a qubit measuring a generic system. The results include the case when a postselection on 
the measured system is made. We provide the exact solution, and also a controlled expansion 
in the coupling parameter, giving formulas valid in the weak measurement regime for arbitrary 
preparation and postselection. The concept of generalized Wigner functions is introduced. 



* dilorcnzo@infis.ufu.br 



1 



I. INTRODUCTION. 



Two-level systems are important because on one hand they realize the simplest non- 
trivial quantum systems, and on the other hand because they can be used as qubits, the 
basic units of quantum computation. While studies on quantum computation have relied 
mostly on projective measurements, it is interesting to consider more general measurements. 
The most general measurements are the positive-operator valued measures that reduce, 
in particular, to postselected weak measurements, introduced in Ref . [7] , for a weak enough 
coupling between system and probe, if the former is subsequently measured projectively. 

Weak measurements are not just a theoretical tool, but they may be useful: for instance, 
it has been shown how weak measurements allow the reconstruction of the wave function [S] 
or of the Bohmian trajectories [9] 

Weak measurements rely on a perturbative expansion of the propagator, and the posts- 
election implies that, according to Bayes' rule, the joint probability of output and postse- 
lection be divided by the marginal probability of successful postselection. While this latter 
probability can be expanded in a Taylor series, its inverse may not, when the preparation 
and the postselection are nearly orthogonal, making the lower order term zero. Hence, one 
should be more careful in making a perturbative expansion in this case [10-13]. This expan- 
sion is sometimes called improperly "non-perturbative" [13], but it actually breaks down for 
larger values of the coupling parameter. A more appropriate term would be perhaps "non- 
polynomial expansion" , as the conditional probability can be written as a rational function 
in the coupling [11]. 

Two-level systems, however, are simple enough to allow for an exact treatment, hence 
there is no need to confine ourselves either to the strong or to the weak measurement 
regime. Indeed, sequential measurements of noncommuting variables have been shown to 
allow quantum state tomography, if one stays in the intermediate coupling regime [14] . In the 
following we shall consider measurements followed by a postselection, but treat the problems 
exactly, while we shall provide approximate formulas valid in the weak measurement regime 
for ease of comparison with the previous results. 

Finally, we shall address two-level systems as "qubits", but we shall use the language 
of spin 1/2, not the language of quantum computation, with an exception in Section III. 
The reason for doing so is that quantum computation relies on a special basis, the so-called 
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computational basis, which conceals the invariance of the formulas. Spin 1/2 systems, on 
the other hand, allow to express significant quantities through geometrical expressions that 
are manifestly invariant, as they rely on the scalar products of vectors. 
The main results of this paper are: 

1. Equation (16), which gives the statistics of a qubit when measured by another qubit. 

2. Equation (28), which gives the statistics of a qubit when measured by a general system. 

3. Equation (39), which gives the statistics of a general system when measured by a 
qubit. 

All these results are exact, apply for any coupling strength, for any preparation of the system 
and of the probe, for any postselection of the system, and for any readout basis of the probe. 

II. GENERAL DESCRIPTION OF A MEASUREMENT 

A measurement proceeds in the following way: (i) prepare the detecting system, which 
we call the probe for brevity, in the ready state pp, uncorrelated to the initial state of the 
system to be measured pi, so that the total state is pp ® Pi; (ii) have the measured system 
and the probe interact through a Hamiltonian i/mt; (iii) read the output of the probe in 
a properly chosen basis E{0), where E{0) is a family of nonnegative operators satisfying 
/ dfi{0)E{0) = 1, with p a Lebesgues-Stieltjes measure accounting for a continuous, dis- 
crete, or mixed output spectrum ['"]. The case most commonly studied is the one where 
E{0) = \Q = 0){Q = 0\, i.e. the readout of the detector is ideal, it is an eigenstate of 
a suitably chosen operator Q. In general, however, E{0) are proportional to mixed state, 
E{0) = pr(0)p(0), with pr(0) being the prior probability of observing an outcome O. Of- 
ten, the prior probability is not normalized, J (ip(O)pr(O) = oo, but this is not a problem: 
what matters are the ratios pr(0)/pr(0') that give the relative probability of observing an 
output O rather than an output O', a priori, i.e. when no measurement is being made. For 
instance, a common case is that all output are equiprobable a priori, so that pr(0) = 1. If 
instead the probe has a different resolution r„ for different outputs 0„, pr(0„) oc r„. For 
ease of comparison with the results of the previous literature, we shall consider the case of 
an ideal output, E{0) = \Q = 0){Q = 0\. 
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In a strong measurement i^int, Pp, and E{0) are chosen in such a way that, if initially 
the system is in the eigenstate \A) of A, the detector ends up in p(0 = f{A)), with f{A) 
a given function of A establishing a correspondence between the output and the state of 
the measured system, and in particular f{A) = XA for a linear measurement. This strong 
measurement is obtained under the following conditions 

1. Choose a bilinear interaction iJint = X6(t)AP, where A is the observable of the system 
we wish to measure and P is an observable of the probe. 

2. Prepare the probe in a state \Q = 0) such that \Qa) = exp{iXAp)\Q = 0) are 
orthogonal for different eigenvalues A of A. 

3. Reading the probe in a basis containing \Qa)- 

Notice that Q and P are not being assumed to be position and momentum of a pointer, 
nor to be canonically conjugated variables, as in the following we shall treat the case of a 
finite-dimensional Hilbert space. 

In general, the one-to-one correspondence between output of the probe and state of the 
system does not hold, either because of intrinsic limitations in the preparation and readout 
operation, or because of intentional deviations from the ideal case with the goal of reducing 
the measurement back-action at the expense of the reliability. In the following, we shall 
drop hypotheses 2 and 3. 

Furthermore, a second measurement on the system can be realized subsequently, allowing 
to select the data of the first measurement depending on the outcome of the second one. This 
procedure is known as postselection. The system can be postselected in a pure or mixed state. 
The postselected state Ef is a positive operator that satisfies Ef < 1 and that may be not 
normalized, but this has no consequences as far as the conditional probability is considered 
Prob{0\postselection) = Prob{0&ipostselection)/Prob{postselection). In particular, the 
case when no postselection is made is obtained for Ef = 1. The postselection in a mixed 
state can be achieved either by making the second measurement a Positive Operator Valued 
measurement [16], or by making the postselection a probabilistic process [11, 17]. Notice that 
if we multiply E{ by a positive constant c, keeping cEf < 1 the normalized postselected state 
does not change. The probability of a successful postselection, however, is proportional to 
Tr E{. Hence, an optimal choice is to choose the maximal c, the one that makes the maximum 
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Symbol 


Meaning 


Restrictions 


P 


Write-in variable of the probe 


None 


Q 


Readout variable of the probe 


None 


Pp 


Initial state of the probe 


Non-negative, trace one operator 


PPIf 


Final state of the probe, given a postselection in 


Non-negative, trace one operator 


T 


Spin operator of a qubit probe 


None 


n 


Initial polarization of a qubit probe 


|n| < 1, pp = (l + n-f)/2 


p 


Orientation of P for a qubit probe 


IpI = 1, P = (l-p-f)/2 


q 


Orientation of the readout of a qubit probe 


|q| = 1 


A 


Variable of the system to be measured 


None 


pi 


Preparation state of the system 


Non- negative, trace one operator 


Ef 


Postselected state of the system 


Positive operator, E{ < I 


Pf 


Normalized postselected state of the system 


Pi = Ef/Tr{Ei) 


a 


Spin operator of a qubit system 


None 


nii 


Preparation polarization of a qubit system 


|mi| <1, pi = (1 + mi •6-)/2 


mf 


Postselection polarization of a qubit system 


|inf| < 1, Pf = (1 + rrif • &)/2 


a 


Polarization of A for a qubit 


|a| = 1, A = a • o" 



TABLE I. The symbols used throughout the paper. 



eigenvalue of E{ equal 1. For instance, let us consider the case where no postselection is 
done, yet a fraction x of the data is discarded by mistake or due to a malfunction of the 
apparatus. Then, Ef — {1 — x)l, the normalized postselected state is pf — 1/N with N the 
dimension of the Hilbert space, the same that would be obtained for x — 0, but in this case 
all of the data would be used. 

III. A QUBIT MEASURING A QUBIT. 

For the sake of exposition, let us take momentarily the Z-axis along the direction of 
the spin component being measured, i.e. A — {1 — a^jl. The shift and the rescaling are 
done for convenience, and are not essential for the results to follow. We choose the readout 
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basis of the detector to be along the Z-axis as well.^ Then, for an ideal measurement, the 
operator P must be a spin component orthogonal to (we indicate with aj the operators 
on the system, and tj the operators on the detector). Without loss of generality, we can put 
P = (1 — Tx)/2. It is readily verified that if the coupling constant is A = vr and the ready 
state is = +1), then the measurement is ideal. Indeed the time-evolution operator is 

(1) 

When the system is in the state \az = +), {(Tz = +|W|o"2 = +) is the identity over the 
detector; when the system is in the state l&z = —), {(Tz = —\U\az = —) = t'x, so that, if the 
detector is initially in the pure state pp = \R){R\, \R) = \tz = +) it ends up in \tz = — ) 
^. In quantum information, the operator U in Eq. (1) is known as a controlled NOT: if we 
identify Tz = + with the logic symbol and = — with 1, then U takes to 1 and vice 
versa for = — (i.e., the control bit is in the "true" state 1), while it does nothing for 

CTz = +■ 

A. Fixed coupling, arbitrary probe preparation 

Along the lines of Ref. [18], we consider a weak measurement, not in the sense that A ^ 1, 
but in the sense that the ready state of the detector is left arbitrary, with A fixed to vr. Hence, 
the detector is considered to be initially prepared in a general state pp = (1 + n • t)/2 with 
|n| < 1. Averages over the initial state will be denoted by a bar. Since the only non-trivial 
averages for a spin-1/2 are the spin components, we shall introduce the average vector r = n. 
Furthermore, to keep equations manifestly invariant, we replace the X-direction of the probe 
qubit with a generic unit vector p, and we write A = sl - & with a a unit vector [notice that 
the quantity appearing in the interaction Hamiltonian is {1 — A)/2, which results in a trivial 
systematic shift and rescaling of the readout]. 

In Ref. [11] we have proved that the statistics of the readout can be written in terms of 
the normal weak values 

a^,„ = Tr(i'"Efi>0, (2) 

^ Notice that the two Z-axes are not necessarily the same. Another way to state this is that we consider 

the eigenstates of the measured quantity A, fix arbitrarily the relative phase, and label them |+), |— ); the 

same is done with the eigenstates of the readout operator Q, which are labelled as \tz = +), \tz = — ). 
^ One can choose as well the opposite ready state jr^ = — ). 



U = exp 
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z-(l - az){l - 
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FIG. 1. Scheme of the measurement of a spin 1/2 using another spin 1/2. 



where A is the observable being measured [a^ in our case). For a spin 1/2, there are 
only three parameters: the positive real numbers ao.o = ^ and ai.i = /3, and the complex 
quantity ao,i = ct- These quantities are already written in invariant form. Their geometric 
expressions are 

1 



UJ 



a 



/3 



Tr(Ef)(l + mi-mf), 
Tr {Ef){m[ + iiif + zrrii x nif) ■ a, 
Tr (Ei) (1 — nii-mf + 2mfaa-mi) 



(3) 
(4) 
(5) 



where the preparation of the system and the (normalized) postselection are pi = (l+mi-CT)/2 
and pf = (1 + mf ■ ct)/2, with |mi| < 1, |mf | < 1. 

We consider the readout to be in an arbitrary basis |q : r). The eigenstates of a ■ ct are 
denoted by |±). The joint probabihty of observing an output r and making a successful 



postselection in Ef is 

V{Ei,T) = (+|Ef|+)(+|pi|+)(q : r|pp|q : r) + (-|^f |-) (-|pi|-) (q : r|p-fppp-f|q : r) 
+ (+|Ef|-)(-|pi|+)(q : r|p-rpp|q : r) + (-|Ef |+)(+|pi|-)(q : r|ppp-r|q : r). 

(6) 

The first line can be interpreted "classically", i.e., as if the measured value is a prepossessed 
property of the system: indeed, if the system has the value +, and this happens with 
probability (+|pi|+), it will be postselected in E{ with probability {+\E{\+), and since in 
this case the evolution of the detector is trivial, the probability of observing the output r is 
(q : r|pp|q : r); if, instead, the system has initially the value 1, then the detector is rotated 
by an angle vr around the X-axis, yielding a probability (q : r|'ra.ppf^.|q : r). The second 
line represents interference effects, does not allow a simple interpretation as the one above, 
and is at the origin of the strange behavior of the weak value. 

The four numbers {(7\E{\a') {a'\pi\a) can be expressed in terms of the four real invariant 
parameters u, (3, Re(a) = a', Im(a) = a", 

(+|Ef|+)(+|pi|+) =l[a; + /3 + 2a'], (7a) 

{-\E,\-){-\p,\-)J-[co + ^-2a'], (7b) 

(+|Ef|-)(-|pi|+) =\[u-fi- 2ta"] , (7c) 

(-|Ef|+)(+|pi|-) =\[iu-P + 2ia"] . (7d) 



The probability of a successful postselection is found by summing Eq. (6) over r 
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V{E,) = \[{l + p-n)u; + (1 - p-n)/3] . (8) 

Next, we consider the conditional probability, obtained by dividing the joint probability 
ofEq. (6) byP(Ef), 

^ Quantum mechanics is a local theory, hence, since the detector and the probe do not interact any longer 
after t — 0, the probability of postselection does not depend on what observable of the detector is measured, 
or on whether an observable is measured at all. This is a natural consequence of the invariance of the 
trace and of the factorization of the time-evolution operator for nonintcracting systems. 
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Since both V{Ef,T) and V{Ef) are homogeneous functions of degree one in am,n, we may 
divide both by u, ehminating one parameter. Accordingly, we define the canonical complex 
weak value Ay^ = a/u and the generalized (real) weak value By^ = fi/uj. For pure preparation 
and postselection Byj = {A^l'^, otherwise^ B^j > 1^^^,^. For brevity, Re(A^) = A'^ and 
Im(y4^„) = A'^. We note that when no postselection is made, A^; = {A)i and B^ = {A'^)i = 1. 

As the outcome is binary, the average value of r gives all the necessary information about 
Q. Multiplying Eq. (6) by r, summing, substituting Eqs. (7), and dividing by Eq. (8), we 
find 

/ \ ^ q-P (1 + P-n) + 2AI (q-n - g-pp-n) - 2A;;(qxn)-p - g^g-p (1 - p-n) 

(1 + p-n) + 5^(1 -p-n) • ^ ^ 

This is a general invariant formula that resumes the statistics of the measurement of a qubit 
with a qubit for fixed coupling A = tt. 



Applications 

In particular, let us consider a readout made in the basis g = z, and an interaction along 
X, as in the first example and corresponding to the coordinate system chosen in Ref. [18]. 
Then 

9 4' ^-1-9 4"^ 

In a strong measurement = 1 , f^. = 7^ = 0. In a weak incoherent measurement, the 
detector is prepared in a mixture of eigenstates of r^, yielding < 1 , f^; = 7^ = 0, so that 
(r) = 2A'yjf^/{l + Byj). In a sense, t\ represents the amplification factor. Following this 
consideration, Ref. [ib] divides the output by = 27^ — 1 and obtains an "amplification", 
even though a two-level system is being used. Indeed, when the initial state of the detector 
is very close to the state \fx = +), say in the XZ plane and making an angle e with the 
X-axis (sine = 27^ — 1 in the notation of Ref. [18]), 

/ \ 2A;„sin£ 

(1 cose) (1 - cos 6)5^ ■ ^ ^ 

In the limit e — )■ 0, (r) / sine — )■ A'^. When the preparation and the postselection are nearly 
orthogonal, this limit can be very large. However, for finite e, the approximate equality 



More precisely [ ], B^j ~ |4.^p whenever the eigenvectors of pf with nonzero eigenvalue, |/ : k), and 
those of Pi, \i : k), satisfy (/ : k\A\i : k') / {f : k\i : k') — constant y k, k' . 
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(T)/sin£ ~ A'^ holds only as far as e'^Byj <^ 1. For nearly orthogonal preparation and 
postselection Byj becomes very large (much larger than A^, recall that Byj > \Ayj\'^) and the 
second-order correction in the denominator can not be neglected. In other words, care must 
be taken to apply the limits £ — )■ and TrjpfPi} = w — )■ in the correct order. 

Interestingly, for £ — > vr, (r)/ sine A'^/B^j, so that preparing the detector in the other 
coherent state If^ — —1) one can retrieve the value of B^. Furthermore, the imaginary part 
of A^ can be inferred by preparing the detector in a state with Ty — 1, — — 0- 

A possible application is the following: Suppose we want to establish whether the initial 
state Pi is pure or mixed. We make a postselection in a pure state pf. By changing the 
preparation of the detector, we retrieve A'^, A'^, and B^,. If By^ > {A'^)"^ + {A'^Y, we 
conclude that the initial state of the system is mixed. 



B. Arbitrary coupling and probe preparation 

Let us consider now a fixed preparation and a varying interaction strength. We shall 
momentarily choose a coordinate system where a = z, p = x, and q = z. We take the initial 
state of the detector to be, in the basis of eigenstates of f^, 

pp = -(l + n-f)= , (13) 

^ 1-wl 

with < w < 1, the probability of the probe being initially in the "up" state, < ac < 



y^w{l — w) the coherence, and 6 a phase. Thus n = 2«;[cos (5)x-|-sin {S)y] + {2w — l)z. The 
interaction is 

H^-h6{t)^{l (14) 

If the detector has a continuous output, the regime of coherent weak measurement is obtained 
for A <^ Hence, one can either increase the coherence k (and with it the spread A = 
2^yw{l — w) > 2k, because of the uncertainty relation) or decrease A. When the detector 
is a system with a finite spectrum, the two procedures are no longer equivalent, since k is 

bounded (in the case of a qubit that we are considering, \k\ < 1/2). 

We find the probability of postselection (we switch back to manifestly invariant expres- 
sions) , 

V{Ei) = ^ [(1 + + S''p-n)uj + 2SC{1 - p-n)a" + -^'(l - p-n)/3] , (15) 
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where C = cos (A/2) and S = sin (A/2). 
The average conditional output is 



(r) = [(1 + C^ + S^p-n) + 2SC{1 - p-n)A': + 5^(1 - p-n)B.^] 
X !^2C\-n + S^q-p{l + p-n) -2SC{qxn)-p 

+ 2A'^ [S' (q-n - q-pp-n) + 5C(qxn)-p] 
+ 2< [5C (q-n - q-p) - ^^(qx n)-p] 

- 5^S'^q-p(l - p-n) \. 



-1 



(16) 



For A = vr, the resuhs of the preceding subsection are recovered. 



In the regime where we can expand up to first order in A, 



(r) ~ q-n+ - [{A^ - l)(qxn)-p - A'^{q-p - q-np-n)] . 



(17) 



Hence, by having q, n, and p mutually orthogonal, the average outcome is proportional to 
A'^ — 1; instead, having q and n orthogonal and q and p parallel, the average output is 
proportional to the imaginary part A'^. 

IV. MEASUREMENT OF A QUBIT WITH A GENERAL PROBE. 

This case admits an exact solution, already provided in Ref. [1'.)] for a simple case (a probe 
in a pure Gaussian state interacting instantaneously with the system). This solution has 
been extended to the case of a finite-duration nondemolition measurement with the probe in 
a mixed Gaussian state [ ], the position and value of the maximum output were estimated 
[20], and the solution was also generalized to an arbitrary initial state of the probe [11]. 
Some limiting cases of the exact solution have also been rediscovered several times [10, 21- 
24]. It has been noticed that the exact solution for a spin 1/2 holds also for any operator 
satisfying = 1 [25, 26]. Furthermore, since fl = (1 — A)/2 is a projection operator, the 
results apply as well to yes-no measurements, which are usually treated separately [27]. 

In the following we provide compact formulas for this latter C3iSG, clS db reference for 
future works on weak measurement and to avoid the publication of further duplicated or 
incremental results. The key observation is that [20] 




(18) 
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Thus, for the instantaneous interaction H = ~XS{t)AP, where P is a variable of the detector, 
the exact joint final state for system and probe is simply 



Ps,p 



cos(AP) +isin(AP)i 



cos(AP) -isin(AP)i 



(19) 



Pi <»Pp 

Upon postselecting the system in the state Ef, the conditional state of the probe becomes 

pp|f =V{Ei)-^ Tts {{Et ® l)ps,p} 

=00 cos(AP)pp cos(AP) — ia* cos(AP)pp sin(AP) + ia sin(AP)pp cos(AP) + (5 sin(AP)pp sin(AP) 

(20) 

The probability of postselection is 



P{Ei) = '^^4^ - Im (a) sin(2AP) + ^— ^cos(2AP), 



(21) 



2 V / V / 2 

where the bar indicates an average over the inital state of the probe. After dividing both 
Eq. (20) and Eq. (21) by u, 

cos(AP)pp cos(AP) - iA^ cos(AP)pp sin(AP) + iA^ sin(AP)pp cos(AP) + By, sin(AP)pp sin(AP) 



pp\i 



[1 + P^)/2 - A'^sin(2AP) + (1 - P^)cos(2AP)/2 



(22) 



This expression can be used to obtain the average of any observable O of the probe, for any 
preparation and postselection, and for any coupling strength. 
For example, the average of a function /(P) is 



Trp{/(P)pp|f} =N~U -A'^/(P)sin(2AP) 



+ 



1 + 



2 /(P) + ^^/(P)cos(2AP)^ 



(23) 



with = V{Ef)/u the denominator in Eq. (22). 



Intermission: generalized Wigner functions 

In order to provide the more general result, we need to introduce the generalized Wigner 
functions. We recall that no assumptions are being made on the variable P, in contrast to 
the literature, where it is assumed that it has a continuous unbounded spectrum. Given an 
additional operator O, with eigenstates \0), we define the generalized Wigner function 

WopiO,P)= [ dfi4p){0\P + p/2){P + p/2\pp\P-p/2){P-p/2\0) (24) 

J Dp 
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where the spectral measure of P is IJ>{P), and we defined as resulting from a change of 
variables P = (Pi + P2) /2, p = Pi - P2 

d/x(Pi)d/x(P2) = [ ^'^^' ^ dfi^Pi - P2). (25) 

This equation defines at the same time the spectral measure of P. We remark that in 
general, the domain of integration of p depends on P, unless P has a continuous unbounded 
spectrum. Furthermore, the generalized Wigner functions are real, as the spectrum of p is 
necessarily symmetric. 

Notice that if P is assumed to be a momentum operator, and O = Q to be its conjugated 
variable, then Wq p is but the ordinary Wigner function W. Furthermore, in this case 
Wqp{Q,P) = Wpq{RQ)., but in general Wo^piO,P) 7^ Wpo{P,0). Notice also that if 
[O, P] = 0, then Wq p{0, P) — pp{P, P)5{P — O). There are some fine points arising when 
the spectrum of either operator is degenerate, which are beyond the scope of the present 
manuscript. 

The generalized averages are defined as 

F{0,P) = J df,{0)df,+{P)F{0,P)Wo,p{0,P) (26) 

with /i(0) the spectral density of O. In particular, if P is a function of O or P only, then 
the generalized average reduces to the ordinary one 

no) = Tr [F{d)pp] , P(P) = Tr [P(P)pp] (27) 

End of the intermission 



Now we are in the position of providing the conditional average of any observable O. The 
average of a function g{0) is 



Trp{s(0)pp|,} =JV-'^ -A'„ [<,(0a) -s(0-a)J -<9(0)sm(2AP) 



where 



6 , , = e±*^^Oe^'"^ 



(29) 



is the operator O in the Heisenberg picture. 
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In particular, if we assume that P has a continuous unbounded spectrum, so that we may 
define its conjugate operator Q, the average of a function g{Q) is 



Trp{^7(Q)pp|f } =N-H-Al[g{Q + X)-g{Q-\)\- A^giQ) sin (2AP) 




l + B^g{Q + \)+g{Q-\) 1 - 



(30) 



where the overhne now indicates an average weighted with the ordinary Wigner function. 

Reference [21] studies the limit By^ ^ A^, and writes Eqs. (21), (23), and (30), in terms 
of Taylor series, without realizing the series could be resummed to give a closed form, as 
was done in Ref. [ii]. 

Finally, we note that, as the results presented here are valid for any coupling strength, we 
could have simplified the formulas by the rescaling P — )■ XP. However, we avoided doing so, 
in order for the reader to be able to compare easily these general results with those obtained 
in the weak coupling limit. 

V. MEASUREMENT OF A GENERAL SYSTEM WITH A QUBIT 

We shall now consider the measurement of an observable A. Such kind of measurement 
was used by Kocsis et al. to determine the Bohmian trajectories of photons [9]. The 
interaction is now taken to be ifmt = —Xh6{t)A{l — p • t)/2. As before, the system is 
prepared in a state pi and postselected in E^. The initial state of the qubit probe is 



A. Approximate results: weak measurement regime 

The weak measurement consists in applying perturbation theory to the joint probability 
V{Ef, t) and to the probability of postselection V{Ef) = J2t "^(-^f > then making a Taylor 
expansion, if applicable, of V{Ef,r)/V{Ef). We expand the time-evolution operator up to 




(31) 
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first-order terms in A. Then, 



P(Ef,r)~TrUEf®|q:r)(q:r|) 



1 + iXA 



1 — p f 



X (Pi ® Pp) 1 - i^A 



1 — p • -f 



uj{l + rq-n) 

+ — p-n + r(q-n — q-p + zqx n-p)] + c.c. 



+ - p-n + ^(q-pp-n - q-p)] 



and 



l-A(l-p-nX + ^A2(l-p.n)5, 



(32) 



(33) 



A fundamental consideration is tliat, wliile tlie time-evolution propagator is expanded to 
first order, the second order term arising when multiplying U and W must be retained, in 
order for the probability to be positive-definite [12]. This prescription is at odds with the 
naive Taylor expansion learnt from calculus textbook, where the approximation applies to 
functions that do not have to obey any constraint, contrary to a probability distribution. 
The conditional average output of the qubit meter is thus 



q-n - A(qxn)-pA^ - A(q-n - q-p)^;; - |A^q-p(l - p n)B 
1 - A(l - p-n)A'; + iA2(l - p-n)B^ 



(34) 



Notice that in section III the measured variable is (1 — a ■ ct)/2, but the normal weak values 
am,n are defined in terms of a ■ &. Thus, in order to recover the results of Sec. Ill, you 
should put A^-^ {1- A^)/2 and {B^ - 2A'^ + l)/4. 

In the linear regime, we recover the results of Ref. [28] 

(r) ~q-n-A(qxn)-pA'^ + A(q-p-q-np-n)A;;,. (35) 

As in the previous section, preparing, interacting, and measuring the qubit meter in orthog- 
onal bases provides A'^, while having the measurement and the interaction bases to coincide 
(or better, not to be mutually unbiased, which suffices) provides A^. 

B. Exact results 



Equation (34) is an approximate expression, while the analogue relations in the previous 
sections are exact. In the following, we shall generalize the approach of Kedem and Vaidman 
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[29] in order to find exact expressions for tlie conditional probability of a readout r of the 
qubit probe in the direction q, given that the postselection of the system in the state Ef 
was successful. Born's rule, applied to the evolved joint state of system and probe, yields 

P(Ef, r) = Tr |(Ef ® |q : r)(q : r|) e^^^^^-P'^^/^ ^ g-ai(i-p.x)/2| 

= '^{<i- t\p : ri){p : Ti\^{l + n- f) |p : r2)(p : Tsjq : r)!^,,^^, (36) 

Tl,T2 

where we introduced twice the identity as J2t Ip • ''")(p • ''"I = 1 defined 

..^ = Trs |Efe^^^(i^^^)/Vie"^'^^'"''^/'} • (37) 

In particular, r+ + = oj, and, for pure preparation and postselection, F^^ factorizes as 
^Ti,T2 = 7ri7r25 with 7 the propagation amplitude ■jr = ("^f | exp [iAA(l — r)/2] \ipi). The 
conditional probability is thus 

E^2(l+™-P)rT,r 

As elsewhere, one parameter can be eliminated in the conditional probability, for instance 
dividing numerator and denominator by u 

^, . Eri,r2(q ■ AV ■■ n)(p : rail (1 + n- f) |p : ri)(p : rslq : r)^,,^^ 

with Cri^r2 = rri,T2/r+,+- For pure preparation and postselection, Cr-^,T2 = ^^0*^, with 
c+ = 1 and c_ = (^/'f | exp [— zAA] lipi) / (ipfltpi) . The latter term was dubbed the modular 
weak value in Ref. [29]. Formally, Eq. (39) is an exact solution. However, it relies on 
calculating terms like {ipil exp [— zAA] lipi). If this term could be evaluated efficiently, then 
we would have an exact solution for any probe, not only the qubit probe being considered 
here. Generally, this is not the the exponential of the operator A is difficult to 

estimate, if we except some simple cases, as, e.g., the case = A considered in Section III, 
or, more importantly perhaps, the case of A being the generator of a transformation. 
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